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Abstract 

We explore a dark energy model with a ghost scalar field in the context of the run- 
away dilaton scenario in low-energy effective string theory. We address the problem 
of vacuum stability by implementing higher-order derivative terms and show that a 
cosmologically viable model of "phantomized" dark energy can be constructed with- 
out violating the stability of quantum fluctuations. We also analytically derive the 
condition under which cosmological scaling solutions exist starting from a general La- 
grangian including the phantom type scalar field. We apply this method to the case 
where the dilaton is coupled to non-relativistic dark matter and find that the system 
tends to become quantum mechanically unstable when a constant coupling is always 
present. Nevertheless, it is possible to obtain a viable cosmological solution in which 
the energy density of the dilaton eventually approaches the present value of dark en- 
ergy provided that the coupling rapidly grows during the transition to the scalar field 
dominated era. 



1 Introduction 



One of the most challenging and intriguing problems of cosmology is undoubtedly that of 
dark energy, the marginally dominant negative pressure component that drives the present 
acceleration of the Universe (for reviews see e.g. ^H]). The identification between dark 
energy and the energy of the vacuum (cosmological constant), which, in some respect, may 
look as a "minimal" choice, rises at least a couple of embarrassing issues. The main concern 
is perhaps to explain why the vacuum energy is so small in particle physics units ("fine- 
tuning" problem). Moreover, beside being generically "small", the cosmological constant 
happens to be exactly of the value required to become dominant at the present epoch. The 
latter mysterious circumstance is sometimes referred to as the "coincidence problem" . 

These puzzles may be better interpreted by assuming that the energy of the vacuum 
is, for some unknown cancellation mechanism, exactly zero and by considering in its place 
a dark energy component with a dynamically variable equation of state. A host of such 
models have been studied, ranging from quintessence K-essence jl], braneworlds jSj, 
tachyons [H], chaplygin gas 12, etc.. Effective scalar fields represent, in this respect, a simple 
and well motivated choice, since they are omnipresent in supersymmetric field theories 
and in string/M theory. For instance, a definite prediction of string theory is that the 
gauge and gravitational couplings are not fixed a priori but, rather, related to the vacuum 
expectation value of a scalar field, the dilaton. More precisely, at the tree level in the string 
loop expansion, the vacuum expectation value of the (four dimensional) dilaton (p is related 
to the gauge coupling acuT and to the Planck mass Mp through e"^ ~ M^/Mp ~ ckgut, 
where Ms = -y/2/a' is the string mass and a' the universal Regge slope parameter of the 
string. 

The very existence of a relation between the couplings and the vacuum expectation 
value (VEV) of a scalar field, although promisingly tasting of "unification", threatens 
with possible violations of the equivalence principle as well as unobserved variations of 
the coupling constants. It is common wisdom [5] to assume that the dilaton and the 
other moduli of the theory are given a mass by some non-perturbative mechanism, in such 
a way that their long range interactions are suppressed and their VEVs frozen at some 
phenomenologically reasonable value. An alternative possibility is that, at the level of 
the effective action, the dilaton decouples from the other fields. To be more explicit, let 
us consider a generic effective string action at lowest order in a' [we use the convention 
( — h ++) for the metric] 



Bg{^)R + B^{<P)g^"' d^^d,^ -^Bf{<P)F^+ ... 



(1) 



where R is the scalar curvature and F is the gauge field. Here the dilaton-dependent loop 
effects as well as the non-perturbative corrections are encoded in the "coupling functions" 
Bi{(j)), and a tilde denotes the quantities as measured in the "string frame", the sigma- 
model metric minimally coupled to the fundamental strings. 

In the non-weak coupling region {e'^ > 1), where the loop effects are important, all 
the (relevant) functions Bi{(j)) may extremize at some value (j)ra- This possibility has been 
investigated in Ref. jH] where it has been shown that the cosmological evolution tends to 
push the dilaton toward the value (pm- This is a phenomenologically "safe" vacuum where 
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the massless dilaton decouples from the matter fields. 

By means of a large-N argument, Veneziano has recently suggested JU] that the effective 
couplings Bi are, for the most part, induced by the quantum corrections of the many moduli 
and gauge bosons of the theory. As a result, these functions of (p reach an extremum at 
infinite bare coupling e'^ ^ oo and follow the general behaviour 

BiicP) = Ci + 0{e-^) (e* » 1). (2) 

In this scenario (of the so called "runaway dilaton" ) , gradually decouples from gravity and 
from the matter fields by evolving towards infinity: 0^ = oo. The moduli of the coefficients 
Ci are proportional, by factors of order one, to the number ~ 10^ of independent degrees 
of freedom which have been integrated over. What about their signs? It is natural to require 
a sensible low energy theory e.g. - in our conventions - that both Cg and Cp\>& positive 
numbers. Moreover, in order for the dilaton </> to behave as a "normal" scalar field in the 
limit of (/> ^ 1, the kinetic coefficient C^^ has to be negative. 

This choice was made in 1121 113| where the late-time cosmology and phenomenology 
of this model with a "well-behaved" scalar field were studied in some detail. In particular, 
in it was shown that a residual coupling of the dilaton to some hidden ( -dark matter?) 
sector of the theory can give rise to a final cosmological attractor with both an accelerated 
expansion and a constant ratio between dark matter and dilatonic energy densities. Such 
a type of "scaling" behaviour jl4l I15j . relaxes the "coincidence problem" by interpreting 
the present energy budget of the Universe (dark matter ~ 1/3, dark energy ~ 2/3 of the 
total) as belonging to a stable attractor configuration. Moreover, these scaling solutions 
can lead to a viable late-time cosmology |16j with the acceleration starting earlier (z > 1) 
than in usual (uncoupled) dark energy models but still consistent with the recent Type la 
supernovae data [T7| ITS]. 

The opposite choice for the kinetic sign in > 0, leads to what in quantum field 

theory is called a "ghost" field or, in a more recent cosmological fashion, a "phantom" 
PH] . When considered as dark energy candidates, ghosts can have an equation of state 
parameter w = p/ p < —1 which is not ruled out (rather, slightly favoured [50]) by current 
observations. Ghosts/phantoms prove to have a viable cosmological behavior which has 
been extensively studied in Refs. jl9l 1211 1^ and also constrained with observations j23j . 
Although intriguing as "classical - cosmological" fields, phantoms are generally plagued by 
severe UV quantum instabilities^, so that the fundamental origin of a w < —1 component 
still represents an interesting challenge for theoreticians. 

Recently, it has been shown that a scalar field with a wrong sign kinetic term does 
not necessarily lead to inconsistencies, provided a suitable structure of higher order kinetic 
terms exists in the effective theory [201 • "^^^ basic mechanism resembles that of a Xip"^- 
theory, where the field acquires a non-zero VEV and the energy of the small fluctuations 
around the new minimum is bounded from below. By the same token, in a theory with 
higher order kinetic terms, the field can "condensate" at a non-constant background value, 

^The energy of a phantom field is not bounded from below and this makes the vacuum unstable against 
the production of ghosts + normal (positive energy) fields |24|. Even when decoupled from the matter fields, 
ghosts couple to gravitons which mediate vacuum decay processes of the type vacuum 2 ghosts + 2'y and 
an unnatural Lorenz invariance breaking - cut off of ~ MeV is required to prevent an overproduction of 
cosmic gamma rays |25|. 
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(pQ 7^ 0, which is perfectly stable at the quantum level. The basic scheme is that of a 
Lagrangian of the type p = —d(p^ + d^p^ /m^ where is a canonically normalized (ghost) 
scalar field and m a mass scale. In this case a background value ipQ can, at the same time, 
be quantum mechanically stable (9(/Jq > and act as a negative pressure component 

(f^Vo — "T-^)- A kinetically driven cosmic acceleration has been proposed both in infla- 
tionary 1^ 1221 and quintessence [1] contexts. See Refs. for recent works about ghost 
condensation. 

In this paper we highlight the runaway dilaton scenario by exploiting the possibility 
that string-loop corrections may result in the "wrong" kinetic sign for the dilaton at the 
effective level [i.e. C^f, > in Eq. ((21)]. The low energy theory can still be consistent if 
the dilaton, by a mechanism similar to that in j^HI) "condensates" in a "safe" vacuum 
configuration. We also show that a residual coupling of the scalar field to dark matter 
can lead to a late-time accelerating Universe with the fractional densities of dark matter 
and dark energy that remain constant in time. In section 2 we sketch the origin of our 
model from the low-energy limit of string theory. In section 3 we address the stability of 
quantum fluctuations for a scalar fleld ip whose Lagrangian is a general function of dip^ 
and of Lp itself. We also generalize the late-time attractors studied in ^1 m and derive 
the general functional form of the Lagrangians that allow this type of scaling solutions. In 
section 4 we study the dynamics of the effective ghost/dilaton for various choices of the 
parameters. It is found that in the absence of the coupling Q between the dilaton and 
the non-relativistic matter one can get a viable attractor solution which asymptotically 
approaches a de Sitter phase without violating the stability of the vacuum. We also show 
that there exist well-behaved cosmological scaling solutions if the coupling Q grows from 
nearly zero to a constant value during the transition to the dilaton-dominant era. We give 
conclusions and discussions in the final section. 

2 A string-inspired model 

Let us consider a general four-dimensional effective low-energy string action [30] 

Bg{<l))R + B^{<l))g^''' d^ct>d,<t) - a'V{(t>) + ... 

+ higher order in a' + Sm[(|>,g^J.v,'^i\ , (3) 

where R is the scalar curvature and (j) is the (four dimensional) dilaton field with potential 
y (</.). Here Sm is the action for the matter fields ^'j, which we suppose that they are 
generally coupled to the dilaton. The dilaton-dependent loop effects as well as the non- 
perturbative corrections are encoded, at the lowest order of approximation in the universal 
Regge slope parameter a', in the "coupling functions" Bi{4>). In the weak coupling region 
of the theory, where the tree-level (four dimensional) string coupling gl = e'^ is much 
smaller than one, the functions Bi{(j)) can be expressed by an expansion of the form 

B,{(t)) = e-^ + CO + cie'*' + . . . (e^ « 1). (4) 

As argued in Ref. il2j, we will assume that in the region with 0^1, the coupling 
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functions Bi have the following behavior 

Bi{(P) = a + Die-^ + 0{e-^^), (5) 

asymptotically approaching a constant value as (p ^ oo. Therefore the dilaton gradually 
decouples from gravity as the field evolves toward the region with (j) ^ 1 from the weakly 
coupled regime. In this scenario, the coefficients Cj and are of order lO'^ (the number of 
independent degrees of freedom which have been integrated over) and unity respectively. 
It is convenient to introduce the conformally related "Einstein" metric jllj 

gf,,. = Cg^Bg{(t))g^^. (6) 

In the above relation the constant factor has been chosen so that the Einstein metric 
approximates the "string" one in the limit — > oo. Note that in the runaway dilaton 
scenario © the 4-dimensional Planck mass Mp is asymptotically related to the string 
scale Ms = ^Jlja! through Mp ~ CgM^, which makes sense if one considers that the 
coefficients Cj are expected to be of order 10^. The action in the Einstein frame reads 



Se= i d'^Xy/^ 



— ^i? — — — V{lp) + higher order terms 



+ Sm[Lp,gf,y,-^i] , (7) 



where we have introduced a canonically defined scalar field <p with the dimensions of a 
mass and the Einstein frame potential V{ip) as follows 



^{B'gV B, 



2 \ Bg J Bn 



ed^\ y(^) = cX'M0)- (8) 



Here a prime denotes a derivative with respect to (j) and e evaluates ±1. When the expres- 
sion in square brackets is positive the dilaton behaves as a normal scalar field (e = +1). 
For example it is known that at tree level, despite the positive sign of the kinetic term in 
the string-frame (S^ = Bg = e~^), we have, in the Einstein- frame, a "proper" scalar field. 

In the runaway scenario that we are considering the value of e crucially depends on 
the sign of C<^, since in the (/> ^ 1 limit the expression in square brackets in Q evaluates 
—C^/Cg + 0{e~'^). It is not difficult to imagine that loop corrections give the "wrong" 
kinetic sign to the dilaton, i.e. > 0. This possibility, in principle, doesn't worry us too 
much since, as sketched in the introduction and analyzed in some detail in section 3, higher 
order kinetic terms may ensure the stability of quantum fluctuations by "condensating" 
|26j the ghost field in a stable vacuum. However, if we are to interpret a ghost condensate 
as dark energy, a "fine-tuning" problem, identical to that of the cosmological constant, 
arises. The higher order kinetic terms are in fact generically suppressed by powers of 
the cut-off (string-) scale and therefore are expected to be largely subdominant today. A 
possible way out, again, is to consider the effects of the loop corrections. In particular, 
the contributions of the lightest fields (neutrinos?) should dominate the higher order 
corrections giving effective terms of the type (dip)'^"' / m^^"'~^\ m being the mass of the 
lightest (coupled) field running into the loops. 

In what follows we will concentrate on an (admittedly, ad-hoc) higher order correction 
to insert in of the type (see also |S2 IH2] ) 

higher order terms = ^(5y?)^e^'^/*^^ , (9) 
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where A and A are numbers of order one and m a mass scale. Exponential corrections of this 
type are naively expected if the theory contains "exponentially light" (m ~ Cm-'r Dme~'^ = 
moe""^; i.e. = 0) fields running into the loops. 

While considering as our basic structure the one with (dip)'^ + (dip)'^, for the sake of 
generality we also cite the possibility of a "zeroth order" -kinetic contribution: a potential 
term V for the dilaton. Since there is no trace of it at tree level, V must be of non- 
perturbative origin (e.g., ~ e"^/^'') and dying out exponentially at weak coupling {<j) — > 
— oo). We will make also the "phenomenological" request that it goes to zero in the limit 
(j) ^ oo. A viable ansatz is 

Vi<P) = Ml [exp(-e-'^//3i) - exp(-e-'^//32)] , (10) 

where Mq is some mass scale, and f3i and [32 are constants satisfying < /32 < This 
is a bell-type potential which has a maximum in an intermediate regime. Note that in the 
strongly coupled region {(j) ^ 1) the potential behaves as 

In summary, the most general scalar field Lagrangian that we consider is 

^dilaton = lid^f + —,{d^)'e'^/'-'- -(w-^) Mle-'^^/^^-, (12) 
z \p2 Pi/ 



where, by equation (jEJ, Ai ~ ^JCgjC^ . Before going into the details of our model, let us 
have a look at some general properties of a scalar field Lp whose Lagrangian density is a 
generic function of {d'-p)^ and of itself. 



3 A (sufficiently) general cosmological scalar field 

In this section we analyze some general features of a canonical cosmological scalar field 93 
whose lagrangian is a general function of the scalar quadratic kinetic term d^"^ and of the 
field itself HZ] . We consider the following action written in the "Einstein frame" 

^Sm\^,^i,g^,y\, (13) 

where X is defined as X = —g^^^d^'pdi,ip/2, and Sm is the action for the matter fields, 
which is generally dependent on ip as well. Hereafter we will use the unit Mp = 1, but we 
restore the Planck mass when it is needed. The Lagrangian (|12)1 can be expressed as 

p(X,v9) = -X + cie^-^X^ -cae-^i'^, (14) 

where ci = AA/m'^ and C2 = (1//32 - l//3i) M^. 

The energy momentum tensor of the scalar field that we obtain from Eq. ()13|) reads 

T^jS) = — 7==T% = 9^luP + PX df^ipduip . (15) 



Sgrav ~\~ 'Sip ~\~ 



R + piX,ip) 
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Here and in the following a suffix X oi (p indicates a partial derivative with respect to X 
and ip respectively. We note that the energy momentum tensor (|15() of the scalar field is 
that of a perfect fluid, T^^ = {p + p)u^Uy + QfiuP-, with velocity = dy^pj \f2X and energy 
density 

p = 2Xvx-V- (16) 



3.1 Consistency conditions 

In order to address the problem of the quantum stability of the scalar field in action (|13() 
we consider the dynamics of the small fluctuations bp around a background value p^ which 
is a solution of the equations of motions. Since p depends both on p and its derivatives, 
the most general case to be considered is that of a background with both (/Jq 7^ and 
Pq 7^ 0. In particular, the choice of a background field evolving in time is essential for 
the mechanism of "ghost condensation" j26j . by which a scalar theory with the "wrong" 
kinetic sign settles to a stable vacuum. 

We divide the field p into a homogeneous part p^ and a fluctuation bp^ as 

p(t,^) = p^it) ^ bp(t,^) . (17) 

Since we are most concerned with UV instabilities it is not restrictive to consider a 
Minkowski background metric and to choose, at least locally, a time direction in such 
a way to satisfy ()17|). By expanding p(X, p) at the second order in bp it is straightforward 
to find the Lagrangian and then the Hamiltonian for the fluctuations. The Hamiltonian 
reads 

n = {px^ 2Xpxx) ^ +Px^-^ - • (18) 

The fluctuations of a time-varying scalar typically obey Lorentz-violating dispersion rela- 
tions This should not worry us too much since, in a cosmological setup, Lorenz in- 
variance is always violated by the presence of a "preferred" (CMB-, comoving observers-) 
frame. The positive definiteness^ of the Hamiltonian is guaranteed if the following condi- 
tions hold 

ii=px + 2Xpxx > 0, ^2 = Px > 0, (19) 

^3 = -Pw > 0- (20) 

Both ^1 and ^2 ai'e usually taken into account when considering the stability of classical 
perturbations. The quantity often used is the speed of sound Cg defined by [SEl 

c?-^ = ^ (21) 

which in cosmological perturbations appears as a coefficient of the k"^ /a? term (here a is 
the scale factor and k is the comoving momentum). While the classical fluctuations may be 
regarded to be stable when is positive, it is essential, for the quantum stability, that both 

^Of course the overall sign of the Hamiltonian for 5if has a meaning if compared with the sign of the 
Hamiltonian of the other fields, which we take to be positive. 
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^1 > and ^2 ^ 0. These two conditions, in fact, prevent an instability of a very bad type 
related to the presence of negative energy "ghost" states which render the vacuum unstable 
under a catastrophic production of ghosts and photons pairs [241 125j . This is essentially a 
UV instability, since the rate of production from the vacuum is simply proportional to the 
phase space integral on all possible final states and since only a UV cut off can prevent the 
creation of modes of arbitrarily high energies. When the Lagrangian p{X, cp) is written as 
p{X, if) = —X — V{(p), as was considered in most of the phantom dark energy models |21j . 
the coefficients and ^2 are both negative. 

The situation is different if we take into account higher-order terms such as X"^ in 
p{X, if). For example, consider the case withp = -X+X'^ 26 . Then one has ^1 = — 1-I-6X 
and ^2 = —1 + 2X. When ^1 > and ^2 ^ 0, corresponding to X > 1/2, the system is 
completely stable at the quantum level. In the region of < X < 1/6 one has ^1 < and 
.^2 < so that the perturbations are classically stable due to the positive sign ofc^. This 
vacua is, however, generally quantum mechanically unstable. 

The kind of instability prevented by the condition H2U() . is of the tachyonic type and 
generally much less dramatic, provided that the conditions (|19|) are satisfied. The point 
is that if pip^p > there are modes, the ones with k'^ < p,p,p/px, that undergo a classical 
exponential growth. This is an intrinsically IR instability of the kind often encountered 
in cosmology when perturbations of super-horizon scales are considered. In fact, for a 
"normal" cosmological scalar field p = X — V{ip), one has just — Pt^^ = V^i^ ~ Hq where 
Hq is the Hubble rate. For the above reasons we will adopt (fT!H) but not (|^n|) as fundamental 
criteria for the consistency of the theory and we will eventually make sure that the effects 
of ^3 are negligible on the relevant physical scales. 



3.2 Scaling solutions 

We now set the field ip in a. cosmological framework and study its evolution in a spatially 
flat Friedmann-Robertson- Walker (FRW) background spacetime with scale factor a{t) 

ds^ = -dt^ + a^{t)d-x^ . (22) 

The equation for 93 can be written as 

V3 {px + ^"^Pxx) + 3Hpx ^ + 2Xpx^ -Pifi = -cr, (23) 

where H = d/a is the Hubble rate and the scalar charge a corresponds to the coupling 
between the matter and the dilaton and it is defined by the relation 6Sm/S(p = —\J—g o. 

We are now interested in a late-time attractor of the kind studied in Refs. jl4l I15L ITT] 
i.e. characterized by a constant $7^ and a constant (and, possibly, positive) acceleration. 
Thus we assume that the Universe is filled by the scalar field (/? and by only one type of 
matter, of energy density pm and general equation of state Wm- In the following we will 
also specialize the formula to the most relevant case Wm = 0, i.e. that of cold dark matter. 
By rewriting Eq. ()23() in terms of the energy density p of the scalar field and then the 
equation for pm we obtain the system 

^ + 3(l + u;> = -Qp„^ 



+ 3(1 + u;^)p = , (24) 



3(1 + Wm)Pm = QPm^ , (25) 
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where 

A^ = lna, w^=p/p, Q{ip) = a/pm- (26) 

We also define in the usual way the fractional density of both components, 0^ = 
Pm/i'iH'^) and Vl^p = p/{3H'^), with + = 1- We are looking for asymptotic scal- 
ing solutions where both the equation of state parameter and the fractional density 
ri(^ are constant during the evolution. In other words this corresponds to dlog/o/dA^ = 
dlogpm/dA^, and by assuming that neither Q is time- varying in the scaling regime we 
obtain, from (|24|) and (|25|) . the relation 



dip 3Qip 



{wm — Wip) = const. (27) 



diV Q 

Inserting Eq. ()27l) back into Eqs. (|24() and H25() . we get the scaling behavior of p and 

dlogp dlogpm on , \ - I o ^ \ ^oQ^ 

diV ^ diV =~3(l + ^^s), with = + 5V(^¥' ~ ^m) • (28) 

Here Wg is the effective equation of state for each component and for the Universe itself 
along a scaling solution. Note in fact that, because of the coupling terms in Eqs. ()24() and 
(|25jl , p and p^ do not scale according to and Wm ■ By the definition of X one also finds 
2X = H'^{dip/dNY oc H"^. Therefore X scales in the same way as p and pm, which means, 
by Eq. that 

^ = -3(l + u;.)X. (29) 

Of course, since the pressure density p = w^pp has a same scaling behavior as p, one has 
dp/dN = —3{1 + Ws)- Since p is also the Lagrangian of the scalar field and it is a function 
of X and if, this implies, by using Eqs. and ((^ . that 

dlogp iSlogp 1 + 0^(7/;^ - -u;^) 

7^1 F~T~7 = A = Q — — ^ ^. (30) 

dlogA A d(p il^{Wm-w,p) 

The solution of this equation gives a constraint on the functional form of p{X, ip) for the 
existence of scaling solutions: 

p{X,ip)=Xg(xe''A , (31) 



where g is any function in terms of y = Xe'^^ . By using Eqs. H27|) and (|29|) one can easily 
show that y is constant along a scaling solution, say, Y = Yq. 

What we have found can be restated by saying that if the Lagrangian of a scalar field 
can be written in the form (|31|) and the scalar is coupled to matter with a coupling Q 
which is weakly dependent of then a scaling solution may take place at late times. The 
"geometrical" properties of such a scaling solution are determined by the parameters Wm, 
A and Q and are independent of the actual functional form of 5. In fact, by the definition 
of A in Eq. I|3U|) . the effective equation of state and the acceleration parameter —q = dajcP' 
are given by 

= A + Q ' -'^ = 2— = 2(ATg^ • 
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Other quantities such as and fi^, that also characterize the late time behavior of our 
model, depend on the actual form of the function g. As already noted, the argument of 
g remains constant along the scaling solution, Y = Yq. A (last) relevant parameter to be 
introduced is therefore the log-log derivative of g{Y) calculated on the attractor. 



d\ogg{Y) 



dlogY 



(33) 

Y=Yo 



Once the functional form of g is given /? can be expressed as a function of the other 
parameters, as we will see in the following with some example. The remaining quantities 
to be determined are then given by 

1 o Q{l+Wm){l + 2P) 

" 1 + 2/3 ' - (A + Q)[wm{l + 2/3) - 1] • ^^^^ 

3.3 Scaling solutions with coupled dark matter (the case Wm = 0) 

Our main interest is the case in which the Universe is filled for the most part with a non- 
relativistic (dark) matter (wm = 0) and the scalar field (p. Since the dilaton follows the 
behaviour of the other component, in the limit Q — > it behaves as a pressureless matter. 
If we assume a non zero interaction between the dilaton and the matter {Q ^ 0), the 
scaling solutions can exhibit the acceleration of the Universe with nonzero 0,^. It may be 
useful to specialize the relevant formulas of last subsection to the Wm = case. Equation 
reduces to 

Q 2Q-X 

"'^ = ^-^- = -aTq' -^ = ^(qTa)- ^ ^ 

While the value of /? depends on the functional form of g{Y), we note that the actual 
numerical value go of g along the attractor is constrained to be 

ffo = -^Q(A + Q). (36) 

The above expression has been obtained by noting that p = SH'^w^Q,^, by using the first 
equation in to express the product WipQip and, finally, by using the relation between 
X and 

;,.,V2^^f^V^^. (37, 



2 \dN J 2(A + g)2' 
The two remaining quantities are obtained directly from Eq. H34|) : 

Since A is positive we get a positive acceleration for Q > A/2. The case Q < —A corresponds 
to a superinflationary attractor (with —q > 1) which we will not consider here. 
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3.4 Simple cases 

It is useful to consider a couple of examples. Let us start with a Lagrangian of the form 
p = f{X) — V{(p). Then Eq. (|5n|) constrains p to be 

p(X, ^)=eX- C2e-^'P , (39) 

where the constant e is either positive or negative. This corresponds to choosing the 
arbitrary function g{Y) as g{Y) = e — C2/Y. With a field redefinition we can always set 
e = ±1 without loss of generality. The scaling solutions for e = 1 (normal scalar field) were 
already extensively studied in the literature (e.g. jlSllllp . while a negative e corresponds to 
the (quantum mechanically unstable) phantom scalar field with an exponential potential. 
The effective equation of state and acceleration of the Universe along the attractor are 
already given by eq. (|35j) once Q and the slope of the potential A are given. In order to 
obtain $7^ and lo^ from l\'68\i we need an expression for the parameter /3. By Eq. we 
have P = (c2/lo)/(e — C2/Yq). The value Yq is related to go by C2/I0 = ^ — go and go is 
constrained by Eq. We finally obtain 



3e + 2Q(A + Q) Q{X + Q) _ + QA + 3e 

2Q{X + Q) ' Q2 + QX + 3e' " (A + Q)2 



When e = 1 these formulas are in agreement with those found in Refs. |151 lllj . Moreover, 
in the non-accelerating case with Q = 0, one has Wip = and = S/X"^, as found e.g. 
in 1^3]. For the phantom type scalar field (e = —1), the fractional energy density fi^ is 
negative for Q = (0^ = —S/X^) which, by itself, means that this case can hardly be 
considered as realistic. This situation may be improved by implementing the coupling Q, 
although we need to caution that the system H39|) with e = — 1 is unstable at the quantum 
level unless higher-order terms in X are taken into account. 

As a second example we consider a higher-order kinetic correction to the usual kinetic 
term, 

p{X,ip)=ciX^e^^-X, (41) 
i.e. g{Y) = ciY — 1. By proceeding as before we obtain 



3-2Q(A + Q) QjX + Q) _ l-QjX + Q) 

2Q{X + Q) ' 3g(A + Q)-3' " {X + QY 



In this case, the stability region is characterized by 2ciXe^'^ > 1 from Eq. (|19() and we 
have a stable late time attractor for X> Q — 3/(4(5). 

3.5 A kind of gauge symmetry 

We conclude this section by making a few comments on a "gauge symmetry" which under- 
lies our equations. The theory of a scalar field with a Lagrangian p{X, ip) as in (|13|) is, in 
fact, invariant under a general redefinition of the field itself. This formally means that all 
the equations have to be invariant under the transformation 




(43) 
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Equivalently, there are different Lagrangians that give the same physics. Take for instance 
p{X, if) and ip) defined as 

p{X,^)=p[x{X,^),^{ip)). (44) 

They must be equivalent since p is obtained from p by a simple redefinition of the field. 
All the evolution equations are in fact invariant under the transformation H43|) . Among the 
positive definiteness conditions for the Hamiltonian in Eqs. and (|2()|) . only (|19|) are 
invariant under H43(l . while H20|) is not. This reinforces the choice of not considering ()2U[) 
as a basic requirement for the quantum stability of the system. 

Under the field redefinition H43|l the scalar charge Q defined in ()26|) varies according to 

QM=Q(^)^. (45) 

dip 

Therefore, when considering scaling solutions in section we "fixed" the gauge freedom 
by taking the charge Q as a constant function of In other words, taking Q independent 
of if is not a restriction, but rather, identifies the scalar field (/? that we are speaking of. The 
functional form of the Lagrangian has thus to be intended as referred to that particular 
"gauge fixed" scalar field. Of course, there is a residual gauge freedom: a rescaling of the 
type 99 — > a with a a positive number. The parameters of the scaling solutions change, 
accordingly, asQ— >Q/a, A— >A/q;, j3 ^ (5 and formulas (|32]) and (|34|) are invariant 
under such a rescaling. 



4 A phase space analysis of various cosmological scenarios 

In this section we shall study the cosmological evolution for the string-inspired Lagrangian 
(fTl)) . The dynamics of our cosmological ghost-dilaton model is best studied by rewriting 
the equations as an autonomous system. We note that, in the special case where Ai = A, 
the Lagrangian (|14() is included as one of the scaling solutions with the choice 

g{Y) = -l + ciY -C2/Y . (46) 

The basic equations we need are the Friedmann equation 

and its first derivative 

2H = ^^- cie^^if^ - pm . (48) 

The latter is equivalent to the evolution equation (j^SJ- 

In order to write the above equations with an autonomous system it is convenient to 
introduce the following quantities: 

x^ = ^ .2 = ^ (49) 
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by which Eq. H47|) can be rewritten as 



with 



2,0 4 / 2 i 2 

X + Scix I z + C2Z . 



(50) 



Fohowing Ref. [23] we shall rewrite the evolution equations in terms of x and z. We assume 
that the matter is non-relativistic (p = 0). Then from Eqs. H47() and (|48() we obtain 



Once the initial conditions of x and z are known, we get the values of and by solving 
Eqs. dSH) and dSSl). 

Basically our main interest is the case ci 7^ and C2 = 0, since the stability condition 
()19() is satisfied even for C2 = 0. This is a minimal string- inspired scenario which ensures 
the stability of the system. In order for the generality of our discussion, however, we 
also implement the C2 term corresponding to the exponential potential of the dilaton. If we 
impose the stability condition (|2flj) in addition to (|19j) . it is required to take into account the 
nonvanishing potential (c2 7^ 0). Although the condition (|2U() is not obligatory compared 
to (fTOj) in order to ensure the stability of quantum fluctuations, we shall discuss the effect 
of the exponential potential in the later subsection. Hereafter we consider the case A = Ai. 

In subsection 14.11 we analyze the parameter range in terms of (2:^,2:^) in which the 
stability of quantum fluctuations is satisfied. In subsection 14.21 it is shown that a viable 
cosmological scenario can be obtained for (5 = without violating the stability of the 
system. This is one of the most important results of our work. In subsection 14.31 we 
shall study the case with nonzero Q and show that it is difficult to get an appropriate 
cosmological evolution unless Q is a rapidly growing function. In Appendix we shall discuss 
the classical dynamics of the phantom field in the absence of the higher-order kinematic 
term (c\ = 0). 

4.1 Restrictions on the configuration space 

In this subsection we derive the conditions that our new variables need to satisfy in order 
to ensure the stability of the system. From the condition ((TO]) we obtain 



where Eq. (|69|) is used. Note that this is also written as y = x^ j z^ > l/(2ci). In what 
follows we shall consider two cases (A) C2 = and (B) C2 7^ separately. 




(51) 



Making use of Eqs. (|^ and (EU, we find 




z^ < 2cix^ 



(54) 
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4.1.1 Caseofc2 = 

When the potential of the dilaton is absent (c2 = 0), the requirement of the condition, 



+ ?>c\x I z < 1, yields 



3cix 

x2 + 1 



(55) 



Note that 0^ > is automatically satisfied under the condition of Eq. H54|) . In Fig. ^ 
we plot the region characterized by Eqs. 1)54(1 and (|55|) in the (x^,z^) plane. The two 
parameters and need be inside this region in order to ensure the stability of system. 
Note that in the case of C2 = the condition ((20} is not satsified. However this condition 
is not fundamentally required for the stability of quantum fluctuations. 

4.1.2 Caseofc2/0 

The constrained region changes in the presence of the dilaton potential (c2 7^ 0). In addition 
to the condition (|54() we obtain the following relation from the requirement 17^ < 1: 



1 

2^ 



12C2X^ 



<z'< 



2c2 L' 



+ 1 + V(a?^+T)2^n^ 



■ (56) 



The stability region ()54j) is plotted together with ()56|) in the right panel of Fig. ^ for C2 = 1. 
Unlike the case of p = X — C2e^^^ in which the constrained regime is + < 1 j33j . the 
allowed phase space is more restricted. Note that there exist allowed regions even in the 
presence of the condition (|2U|) . i.e. \J cxjc^ x^ < z^. This happens for nonzero values of C2, 
but Eq. (|^n)) is not obligatory as we already mentioned. 



4.1.3 Equation of state 

We wish to point out general features when we impose the stability condition l|19l) . The 
Px + ^'^Pxx term in Eq. (|23|) is positive as long as Eq. (|19|) holds. This means that 
the dilaton does not climb up the potential hill unlike the case of the classical evolution 
discussed in the Appendix. The equation of state for the field f is given as 



w,. 



^ X dp 
'^dX 



1 



(57) 



In order to achieve the equation of state with = p/ p < 0, we require the negative 
pressure {p < 0). If we impose the condition (|19|) . one has 2{X/p){dp/dX) < and 
> —1. Therefore as long as the stability condition (|19() is imposed, one gets the 
equation of state for the normal scalar field, in spite of the fact that the coefficient of the 
ip'^ term is negative. This is a general feature of any form of the Lagrangian p. 



4.2 Dynamics of the system for Ci 7^ and Q = 

Let us consider the dynamical evolution for the system characterized by the Lagrangian 
(fTl)) with (5 = 0. In the next subsection we shall study the case in which the coupling Q 
is taken into account. 
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Figure 1: The region in which the conditions H19|) and < 1 are satisfied for (A) ci = 1, 
C2 = and (B) ci = 1, C2 = 1. In the case (A) each curve or hne corresponds to (i) 
z'^ = (3cix^)/(a;^ + 1) and (ii) = 2cix^, whereas in the case (B) each corresponds to (i) 

v2 



1 

2C2 



X 



+ 1 ± y/{x'^ + 1)2 - 12c2x4 , (ii) = 2cix'^ and (iii) = y^ci/c2 
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The discussion in Sec. 13. 21 demonstrates that one has a scahng solution where p is 
vanishing for Q = 0. This is not an interesting attr actor solution, since it does not lead 
to an acceleration of the Universe. As studied in Appendix, it is expected that there exist 
other attractor solutions where approaches 1. These solutions can be found by setting 
dx/dN = and dy/dN = in Eqs. and Especially we are interested in the case 
where w^p approaches a negative constant close to —1. From Eq. (|57|) the equation of state 
for the Lagrangian p4j) is given by 



-l + 3ciY + C2/Y 



Then one may look for the attractor solutions which keep the quantity Y 
Setting dx/dN = and dy/dN = in Eqs. (O and (|^, we get 



VG{2ciY - 1) 



A(3ciy - 1 + C2/Y) ' 
Eliminating x from these equations gives 

24c?y^-3c?(16 + A2)y4 



and x'^ 



1 - \/6Ax/3 
I-C1Y + C2/Y 



x^ / constant. 



(59) 



(24c? C2 - 30ci 



+(24ciC2 - 6 - A^ + 2A2ciC2)y2 



6c2y + }?(?2 



4A2ci)y3 
0. 



In the limit of A — > one obtains the following solutions 



y = 0, 



1 1 + VI + 4ciC2 



(60) 



(61) 



2ci ' 2ci 

where we used Y = Xe^"^ > 0. The ghost condensate considered in Ref. P2] corresponds 
to the attractor Y = l/(2ci) with C2 = 0. Generally we have two stable solutions, i.e. 
y = l/(2ci) and y = (1 + vT+~4ciC2)/(2ci). The former one is an attractor solution 
in which the system evolves toward the de-sitter phase characterized by a cosmological 
constant {w^ = —1). The latter one is the scaling solution for Q = discussed in Sec. 13.21 
but this does not correspond to an accelerating universe since the equation of state is 



w,, 



0. 



In the case of nonzero A, it is generally not easy to get an explicit form of solutions for 
Eq. (|6fl)) . However we can obtain numerical values of Y once ci, C2 and A are specified in 
Eq. (jnOJ. In what follows we shall consider the cases (A) C2 = and (B) C2 / separately. 



4.2.1 Case of C2 = 

When C2 = Eq. (|^ gives 
1 1 

Cl ' Cl 



y = 0, 




1 

Cl 




(62) 



Among them the solutions satisfying the condition ()19() are the second and third ones. 
The former one is a non-accelerating scaling solution with Wip = 0. In the latter case the 
equation of state (|58|) is written in terms of A: 



1 - (AV8) - (A/24)V9A2 + 48 
1 + (3A2/8) + (A/8)V9A2 + 48 



(63) 
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Figure 2: The variation of i^m, Q^p, w^p, and for ci = 1, C2 = and A = 0.1 with 
initial conditions Xi = 0.0085 and Zi = 0.0085. The solution evolves toward the attractor 
characterized by constant and z"^ with — > 1 and — ^ 0. The equation of state 
asymptotically approaches a constant value Wip ~ —0.889. 



In the limit of A ^ we have w^p ^ — 1. The equation of state ranges in the acceleration 
region — 1 < w^, < —1/3 as long as A < 0.82. For example one has = —0.889 for 
A = 0.1. 

In Fig. 121 we plot the evolution of ^ip, w^, and z^ for ci = 1, C2 = and 
A = 0.1. The solution evolves toward the attractor characterized by constant x^ and z^ 
with $7^ — > 1 and — *■ 0. The third solution in Eq. (|62|) gives Y ~ 0.515 for ci = 1 
and A = 0.1. Making use of Eq. and Y = x^ j z^ ~ 0.515, we obtain x^ ~ 1.834 and 
z^ = 3.561, in which case il^ = —x^ + 3cix^/z^ asymptotically approaches 1. This is 
different from the scaling solution that approaches a nonzero value of Q.^^ or Qrn discussed 
in Sec. l3.2l Note that the discussion in Sec. 13. 21 does not apply to the case in which Vim or 

approaches zero asymptotically. 

Numerical integrations of Eqs. H52|) and (|53() agree very well with our analytic estimates, 
see Fig. 121 We begin to integrate the equations with initial conditions satisfying x^, ^ 1 
from the epoch where matter dominates and radiation can be neglected. The case shown 
in Fig. [21 corresponds to the initial matter density — 9.0 x 10~^^^Mp, in which case 
grows to 0.7 at the present epoch ~ 4. The system asymptotically evolves toward the 
values = 1, Q.^ = 0, x^ = 1.834 and z^ = 3.561. We found that the results hardly 
change even if the radiation is included at the beginning. We also checked that the system 
always lies inside the quantum mechanically allowed region shown in the left panel of Fig. ^ 
This is also confirmed in Fig. (21 in which the equation of state is larger than —1. 
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Figure 3: The variation of Jim, f^</3, '^^(/j! and for ci = 1, C2 = 1 and A = 0.1 with 
initial conditions Xi = 0.003 and Zj = 0.0035. The evolution of and $7^. is similar to 
Fig. 121 although the asymptotic values of x^, and Wip are different. 



4.2.2 Case of C2 / 

Let us next consider the situation where the exponential potential of the dilaton is present 
(c2 7^ 0). Equation (|60j) is more involved than in the case of C2 = 0, but we can numerically 
solve Eq. (|?Tn|) for given values of A, ci and C2- For example, when A = 0.1 and ci = C2 = 1, 
the solutions are Y ~ 0.00167,0.467,0.532, 1.62. Since the last one {Y ~ 1.62) corresponds 
to the non-accelerating universe, only the third one {Y ~ 0.532) is an important solution 
satisfying the stability condition (fTU]) . In this case the equation of state is given as ~ 
—0.948 by Eq. (|58j) . which is slightly larger than —1. We also obtain the asymptotic values 
x'^ = 0.4039 and z^ = 0.7591, in which case the fractions of the kinematic energy and the 
potential energy of the dilaton are Qk = — x^+3cix^/z^ ~ 0.2409 and = C2^^ — 0.7591, 
respectively. Therefore = i^x + asymptotically approaches unity. The numerical 
simulation in Fig. |31 agrees well with these asymptotic values of w^p, x'^, and Q.^p. 

By taking into account the nonzero C2 term, the allowed region in the (x^,z^) plane 
tends to shift toward smaller values of x^ and as seen in Fig. Q This does not necessarily 
mean that the cosmological evolution is more restricted than in the case of C2 = 0. In the 
case shown in Fig. |31the system always lies inside the quantum mechanically allowed region 
shown in Fig. ^ In addition even the condition (|2()j) is satisfied throughout the evolution. 

In summary we can obtain stable cosmological solutions which asymptotically approach 
r^ip = 1 and f^m = in the future both in the case of C2 = and C2 / 0. We stress that 
this is different from the scaling solutions discussed in Sec. 3. 
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4.3 Dynamics of the system for Ci 7^ and Q ^ 



If there is a coupling between the dilaton and the matter {Q ^ 0), one has additional 
scaling solutions discussed in Sec. 13.21 When Wm = we have w^, = —Q/{^lip{X + Q)) from 
Eq. (|35j) in the scaling regime. 

Using the relations = p/p, = p^/{3H^), 3H^ = 2(A + QfX/3 and p = X{-1 + 
ciY — C2/Y), Eq. gives the constant value of Y: 



Y 



1 
2^1 



-Q{X + Q) + 



+ 4ciC2 



(64) 



Note that the third solution in Eq. (|61j) is recovered in the limit of Q — > 0. Substituting 
Eq. (jni) for % = p^/{3H'^) with p^ = X{-1 + 3ciY + C2/Y), we obtain 



(A + Q)2 



1-Q(A + Q) + 



12ciC2 



3 - 2Q(A + Q) + V[3-2Q(A + Q)]2 + 36ciC2 



• (65) 



The equation of state is determined from Eqs. (|35|l and (|65|) . 

If one requires the late time behavior with ~ 0.7 and w,^ ~ —1, Eq. gives the 
relation A = (3/7) Q. Substituting this relation for Eq. (|65|) . we can determine the values 
of Q and A for known values of ci and C2- For example, in the case of ci = 1 and C2 = 0, 
one has Q = 0.678 and A = 0.398 for ~ 0.7 and Wip ~ —0.9. As another example, when 



ci = 1 and C2 = 1, Q ~ 1.45 and A ~ 0.85 for Q^:, ~ 0.7 and Wu 



-0.9. 



4.3.1 Case of constant Q 

Let us consider the situation in which the coupling Q is constant not only in the scaling 
regime but also in the matter dominated era. We first analyze the case where the ex- 
ponential potential is absent (c2 = 0). In Fig. |l]we plot the phase-space trajectories for 
Q = 0.678, A = 0.398 and ci = 1 with different initial conditions of and z^. One has 
an attractor point (x^,Zg) = (1.295,2.522) corresponding to the scaling solution given in 
Eqs. (|35|) and (|65|) . In order to ensure the stability of quantum fluctuations, the trajec- 
tories are required to be inside the region shown in Fig. ^ As seen from Fig. |SJ only the 
initial conditions which are not far from the attractor point (x^, Zg) satisfy this constraint. 
There are some trajectories which finally approach the attractor point but are away from 
the stability region in the middle of the way. In this case the vacuum is unstable at the 
quantum level during some moment of time, so we can not regard this as ideal scaling 
solutions. 

During the matter dominant phase and z'^ are needed to be much smaller than unity. 
This means that the trajectories start from the region with a;^ , <C 1 as a realistic dark 
energy model. However, these trajectories do not satisfy the stability condition as seen in 
Fig. |1I In addition they fall into the singularity on the line z'^ = Gcix'^ at which the speed 
of sound diverges. It is inevitable to avoid this divergent behavior as long as x'^ is initially 
less than of order 0.1. 

This situation does not change much even when the exponential potential of the dilaton 
is present (c2 7^ 0). In Fig. [51 the phase-space trajectories are shown for Q = 0.678, 
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Figure 4: The phase-space trajectories for Q = 0.678 and A = 0.398 for ci = 1 and 
C2 = 0. In this case there exist scahng solutions that asymptoticahy approach ~ 0.7 
and ~ —0.9, as long as the initial values of and are not far from the attr actor 
point {x1,z1) = (1.295,2.522) (this point is indicated by a black point in the figure). If 
the initial values of and are much smaller than 1, the trajectories tend to evolve out 
of the stable region. In particular it happens that the solutions hit the singularities with 
= at which the speed of sound diverges. 



19 



0.1 0.2 0.3 0.4 

Figure 5: The phase-space trajectories for Q = 1.45 and A = 0.85 for ci = 1 and C2 = 1. In 
this case scaling solutions correspond to x1 = 0.283 and = 0.505 with asymptotic values 
VL^p ~ 0.7 and w^p ~ —0.9. As is similar to Fig. ^ the trajectories tend to be away from the 
stable region when the initial values of and are much smaller than 1. 

A = 0.398, ci = 1 and C2 = 1. In this case the attractor point corresponds to {x1,zl) = 
(0.283,0.505), whose values are smaller than in the case of C2 = 0. However the behavior 
of solutions with initial conditions x?,^^ <C 1 is similar to the one discussed above, which 
means that the trajectories hit the singularities at = 6cix^ in addition to the violation 
of the stability condition. Therefore inclusion of the constant coupling Q does not help to 
obtain realistic scaling solutions where V^^p is initially much smaller than 1 and evolves to 
a constant value ~ 0.7. 

4.3.2 Case of varying Q 

While the above discussion corresponds to a constant value of Q from the matter-dominant 
to the dilaton-dominant era, one can generalize the analysis to the case in which Q is a 
time-varying function. The reason why we did not get viable cosmological solutions in the 
previous subsection is that the coupling between the matter and the dilaton is too large in 
the matter-dominant era. If Q varies in a way such that it is small in the matter-dominant 
era but grows as the system enters the dilaton-dominant phase, it is expected that we may 
get viable cosmological scaling solutions at late times. The authors in Ref. TT considered 
the coupling of the form 

where Qo ^-nd h are constants. This is suppressed in the weakly coupled regime, but 
approaches a constant value Qo as the field evolves toward large values. 
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Figure 6: The variation of fi^ in the case of the varying Q for A = 0.85 and ci = C2 = 1 with 
initial conditions Xi = 0.003, Zi = 0.0035 and ipi = 6.6. The cases (i) and (ii) correspond 
to the choices {i) b = 7 x 10^, Qq = 1.45 and (ii) 6 = 2 x 10^ Qo = 1-45 in Eq. In the 
case (iii) we choose the couphng Q which suddenly changes from Q = to Q = Qo around 
= 0.7. Inset: The evolution of the equation of state for the cases (i) and (ii). 



In Fig. 121 we plot the evolution of for several different cases with A = 0.85, ci = 1 
and C2 = 1. In the case (i), corresponding to 6 = 7 x 10^ and Qq = 1.45 in Eq. the 
coupling Q is negligibly small during the matter-dominant era. Therefore the dynamics 
of the system is similar to the Q = case before 0^ reaches close to unity. The effect 
of the coupling Q becomes important after that, which eventually leads to the approach 
to the scaling solutions discussed in Sec. 13.21 As seen in Fig. El we numerically obtained 
attractor values $7^ = 0.7 and = —0.9 for the choice A = 0.85 and Qq = 1.45. Note 
that the equation of state always ranges in the region > —1, which means that the 
vacuum is stable at the quantum level. This case corresponds to the solution in which i}^ 
is growing now (N = 4) but approaches a constant value with < $7^ < 1 in the future. 
The asymptotic value of il^ is dependent on the choice of A and Qq. 

The initial values of Q get larger if one chooses smaller values of b in Eq. H66|) . Although 
this is required to make the time before the system approaches the scaling solution shorter, 
the larger initial values of Q lead to another problem. As seen in the inset of Fig. |H1 the 
equation of state w^p becomes less than —1 for the case (ii), marking the instability of the 
system. In order to avoid this problem, the initial value of Q needs to be sufficiently small 
relative to unity as in the case (i) of Fig. El If Q jumps rapidly from nearly zero to the 
value Qq around = 0.7, it is possible to connect to the scaling solutions with ~ 0.7 
around = 4. This is shown in the case (iii) of Fig.O in which we used the step-function 
for the coupling Q. 

If one adopts the function (|66|) . it is not easy to get the scaling solutions which directly 
approach = 0.7 as in the case (iii) of Fig. El The situation is similar when the potential 
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of the dilaton is absent (c2 = 0). This is mainly because inclusion of the larger coupling 
Q tends to work to violate the stability condition H19|) during the transition to the scalar 
field dominant phase. Nevertheless it is intriguing that we found the scaling solutions as 
in the type (i) that evolve to a constant in the future without violating the stability 
condition. 

5 Conclusions 

In this paper we have studied a new type of dark energy model based on string theory. Our 
starting point is the Lagrangian © in low-energy effective string theory. We adopted the 
runaway dilaton scenario in which the dilaton field is effectively decoupled from gravity 
as the field evolves toward i;^ — > oo. The coefficient of the kinematic term of the dilaton can 
be negative in the Einstein frame, which means that the dilaton behaves as a phantom-type 
scalar field. We implemented higher-order derivative terms for the field which ensure the 
stability of the system even when the coefficient of (j)^ is negative. 

We then considered a general Lagrangian which is generic function of a (canonically 
defined) scalar field and X = —{\/ip)'^/2. We first derived the stability condition of 
quantum fiuctuations, see Eq. H19|) . At the classical level the perturbations may be stable 
if the speed of sound given in Eq. (|2H) is positive, but we obtained a more stringent 
condition in order to ensure the stability of the vacuum at the quantum level. When the 
condition (|19|) is satisfied, one can avoid the gravitational creation of ghost and photon 
pairs which was recently raised as a serious problem by Cline et al. j25j for the classical 
phantom models considered so far ^1 1^ [221 [22] . In this work we carefully addressed the 
way to overcome this problem. The stability region is plotted in the phase-space, which 
restricts the dynamical evolution of the system relative to the normal scalar field with a 
positive kinematic term. In Sec. 13. 21 we derived a general form of the Lagrangian for 
a single field ip when the scaling solutions exist. We implemented a coupling between the 
dilaton and the matter (denoted by Q) and obtained the equation of state and the 
relative energy density for the field ip, see Eq. This is a general useful expression 

valid not only for the normal scalar field but also for the ghost-type scalar field. We also 
found that when the Lagrangian is written as p = f{X) — V^(</>) the requirement of the 
scaling solutions uniquely determines the Lagrangian as p = ±X - C2e-^'t'. This is the 
system with an exponential potential and a positive/negative kinematic term, which was 
extensively studied in the literature. 

We considered the classical dynamics of the system with a negative kinematic term 
—X and the bell- type potential (|T?1|) . and reproduced the result that the dilaton evolves 
toward the potential maximum with an equation of state w^p < —1. However, since this 
system is unstable at the quantum level, we took into account a term of the type e^'^X'^ 
in order to ensure the stability of quantum fiuctuations. When the condition for the 
stability is satisfied, the equation of state is restricted in the range Wip > — 1 in general. 
We found that there exist attractor solutions with Q = in which fi^ evolves toward 1 
and Wip approaches a constant value slightly larger than —1. This is different from the 
scaling solutions discussed in Sec. l3.2l which keep the relative energy density constant with 
nonzero values of and Q^n- This can be derived by setting dx/dN = dy/dN = in 
the evolution equations (|52|) and (|53|) . We plotted such examples in Figs. Hand El whose 
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results agree well with analytic estimates. We also checked that the system lies in the 
quantum stability region shown in Fig. ^ throughout the evolution, thereby providing a 
viable dark energy model even when the coefficient of the (fP' term is negative. When the 
matter Lagrangian Sm is not dependent on the field cj) in the string frame, one can show 
that Q oc {1 / B g)(iB g / is vanishing as ^ oo for the choice ©. In this case one can 
get ideal cosmological solutions as shown in Figs. El ^-nd |21 If the matter is coupled to 
the dilaton in the string frame, the coupling Q is not necessarily zero. We analytically 
found in Sec. 13.21 that scaling solutions exist when Q is a nonzero constant in the scaling 
regime. When Q is a constant from the matter-dominated to the dilaton-dominated era, 
the stability condition is violated for the realistic initial conditions {xj, zf <C 1). Therefore 
it is difficult to obtain a viable cosmological evolution in such a case, see Figs. 0] and |SJ 
However, if the coupling Q rapidly grows from nearly zero to a constant value as in Eq. (|(if))) . 
the system can approach the scaling solution without breaking the stability condition, see 
the case (i) of Fig. El In order to get a "direct" approach toward $7,^ ~ 0.7 as in the case 
(iii) of Fig. El we require a step-like change of the coupling Q during the transition to the 
scalar-field dominant era. 

In our work we carried out the analysis for the simplified Lagrangian (|14() in order to 
understand the basic picture of the system with a phantom-type scalar field. In string 
theory we have other non-perturbative and loop corrections such as the Gauss-Bonnet 
curvature invariant |311 1321 136j . It is certainly of interest to extend our analysis to such a 
direction. 
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APPENDIX : Classical dynamics of the phantom field 

In this Appendix we consider the classical evolution of the system in the absence of the 
higher-order terms in X with the exponential potential of the dilaton (ci = 0, C2 7^ and 
Q = 0). Although the stability of quantum fiuctuations is not ensured for ci = 0, it is 
worth clarifying the classical evolution to compare it with the case in which higher-order 
terms are included. When Q = we obtain the following fixed points by setting dx/dN = 
and dz/dN = in Eqs. (l52|l a.nd 
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(67) 
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The first solution is a trivial one which approaches — > at late times. On the other 
hand the second one corresponds to an asymptotic solution with 0,^ ^ 1 at late times. 
This differs from the scaling solutions that approach a nonzero value of or The 
scaling solutions which keep J7<^ constant in the range < < 1 do not exist in the 
present case, since fi^ becomes negative in Eq. (IIU)) for Q = and e = — 1. 
When ci = the equation of state for the field is given as 

= . (68) 

From Eq. ()49|) the function Y = Xe^'^ is generally expressed as 

Y = x^jz^ . (69) 

Therefore the second fixed point in Eq. (|67() corresponds to the equation of state 

w^ = -l-X^/3, (70) 

which is less than —1. It is clear that this is different from the scaling solution H40|) with 
Q = 0, since — > in this case. This constant behavior of for the exponential 
potential was numerically found in Ref. but we succeeded to derive it analytically. 
The field ip is driven up the potential hill due to the negative kinematic term. 

We can consider the bell-type potential in order to follow the classical evolution 
of the dilaton. Then the dilaton evolves toward the potential maximum and is eventually 
stabilized there. In the region where the potential is effectively described by an 
exponential one, we have numerically checked the existence of the constant equation of 
state characterized by Eq. (f7fl|) . For the bell- type potential, this stage is only transient and 

approaches —1 as the field evolves toward the potential maximum. After the field is 
stabilized there, we have x = and 1" = 0, thus yielding = —1. This classical evolution 
of the dilaton is similar to what was already discussed in Refs. |241 123j . hence we do not 
repeat it here. If one takes into account the coupling between the dilaton and the matter 
(Q / 0), we have additional scaling solutions characterized by Eq. (^OJ. 
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